Single top quark production at the Tevatron: threshold 
resummation and finite-order soft gluon corrections 
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Abstract 



I present a calculation of threshold soft-gluon corrections to single top quark pro- 
duction in pp collisions via all partonic processes in the t and s channels and via as- 
sociated top quark and W boson production. The soft-gluon corrections are formally 
resummed to all orders, and finite-order expansions of the resummed cross section are 
calculated through next-to-next-to-next-to-leading order (NNNLO) at next-to-leading 
logarithmic (NLL) accuracy. Numerical results for single top quark production at the 
' Tevatron are presented, including the dependence of the cross sections on the top quark 

' mass and on the factorization and renormalization scales. The threshold corrections in 

OO ■ the t channel are small while in the s channel they are large and dominant. Associated 

I tW production remains relatively minor due to the small leading-order cross section 

^ ' even though the K factors are large. 

^ ' 

o . 

i!B ■ 1 Introduction 

PhI The top quark was discovered in 1995 by the CDF and DO collaborations at the Tevatron 
in top-antitop quark pair production events [H [2]. The search for single top quark events 
continues intensively at the Tevatron [3lll]. The theoretical cross section is less than that of 
the ti cross section and the backgrounds make the extraction of the signal challenging. 
Single top quark production provides opportunities for the study of the electroweak prop- 
5^ I erties of the top quark (such as a direct measurement of the Vtt CKM matrix element), for 
further insights into electroweak theory since the top quark mass is of the same order of 
magnitude as the electroweak symmetry breaking scale, as well as for the possible discovery 
of new physics (extra quarks or gauge bosons, modified top quark interactions, etc.) (see, 
for example, Refs. [5-11]). 

The production of single top quarks can proceed through three distinct partonic processes. 
One of them is the t-channel process that proceeds via the exchange of a space-like W boson 
(Fig. 1), a second is the s-channel process that proceeds via the exchange of a time-like 
W boson (Fig. 2), and a third is associated tW production (Fig. 3). In the t channel we 
have processes of the form qb — > q't and qb — > q't, such as ub dt and db — > ut as well as 
processes involving the charm quark and Cabibbo-supressed contributions. In the s channel 
we have processes of the form qq' bt, such as ud bt as well as processes involving the 
charm quark and Cabibbo-supressed contributions. Associated production of a top quark 
and a W boson proceeds via bg —>■ tW~ as well as Cabibbo-supressed contributions. 

At the Tevatron the t-channel process is numerically dominant, the s-channel process 
is smaller, and associated tW production is quite minor. Although we will be discussing 
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the production of a top quark in this paper we note that the corresponding cross sections 
for the production of an antitop quark at the Tevatron are identical Calculations of NLO 
corrections to single top production via the various partonic modes have appeared in Refs. 
[12-20]. 

The cross section for single top quark production in proton-antiproton collisions is 

^ = XI / dxidx2 0/i/p(xi, hf) (pf2/p{x2, IJ'f) d-{s, t, u, fxp, /iR, Os) , (1-1) 
/ 

where (0/2/p) is the distribution function for parton /i (/2) carrying momentum fraction 
Xi (X2) of the proton (antiproton), is the factorization scale, and is the renormalization 
scale. The parton-level cross section is denoted by a and can be written as a perturbative 
expansion in the strong coupling a^, and s, t, u are standard kinematical invariants formed 
from the momenta of the particles in the hard scattering. 

Near kinematical threshold for the production of a specified final state, such as a single top 
quark, large corrections appear from soft-gluon emission [211 ES, [231 Elj- These corrections 
arise from incomplete cancellations of infrared divergences between virtual diagrams and 
real diagrams with soft (i.e. low-energy) gluons. In principle one can obtain the form of 
these soft radiative corrections at any order in and formally resum them to all orders. 
However in practice such resummed cross sections depend on a prescription to avoid the 
infrared singularity and ambiguities from prescription dependence can actually be larger than 
contributions from high-order terms [24j. Hence, we provide fixed-order expansions of the 
resummed cross section, as has been done for many other processes [TO l [2^ [25l[26l [27 1 l28l[29] . 
in order to avoid such ambiguities. We will calculate soft-gluon corrections for single top 
quark production through NNNLO at NLL accuracy. 
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Figure 1: Leading-order t-channel diagram for single top quark production. 

In Section 2 we briefly describe the threshold resummation formalism and provide expres- 
sions for the resummed cross section. In Section 3 we expand the resummed cross section in 
powers of as and provide formulas for the soft-gluon corrections through NNNLO. In Section 
4 we present numerical results for single top quark production via the t channel at the Teva- 
tron. Analogous results are provided for the s channel in Section 5 and for associated tW 
production in Section 6. The conclusion is in Section 7, and the Appendix collects formulas 
and details on the kinematics and on electroweak parameters used in the calculations. 
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Figure 2: Leading-order s-channel diagram for single top quark production. 



b 



t 



> 



b 




9 



W 



Figure 3: Leading-order associated tW production diagram for single top quark production. 

2 Threshold resummation 

In this section we present the analytical form of the resummed cross section for single top 
quark production. Details of the general resummation formalism for hard-scattering cross 
sections have been presented elsewhere [2T1 [221 [221 [211 [2S1 [30] so here we explicitly show only 
the expressions directly relevant to single top quark production. 

For the process pi + P2 —>■ P3 + Pa (where the particles are represented by their momenta 
Pi), the partonic kinematical invariants are s = {pi +^2)^, t = (pi —Ps)'^, u = {p2 —paY, 
S4 = s + t + u — m\ — m\. Note that we ignore the mass of the 6-quark in the kinematics. 
Near threshold, i.e. when we have just enough partonic energy to produce the final state, 
S4 ^ 0. The threshold corrections then take the form of logarithmic plus distributions, 
[In'(s4/m^)/s4] + , where / < 2n — 1 for the n-th order QCD corrections and rrit is the top 
quark mass. These plus distributions are defined through their integral with any smooth 
function, such as parton distributions, giving a finite result [2T| [30]. 

The resummation of threshold logarithms is carried out in moment space and it follows 
from the refactorization of the cross section into hard, soft, and jet functions that describe, 
respectively, the hard scattering, noncoUinear soft gluon emission, and collinear gluon emis- 
sion from the partons in the scattering [211 [221 [23] . We define moments of the partonic cross 
section by (y{N) = J^ds^/s) e~^'^*^^a{s4), with N the moment variable. Under moments log- 
arithms of S4 transform into logarithms of A^, which exponentiate. The resummed partonic 
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cross section in moment space [2T| [22l [23| [30] is then given by 
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The first exponent resums colhnear and soft gluon emission from the incoming partons 
in the hard scattering and is given in the MS scheme by 
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TT 



0{al). (2.2) 



Here Ni = N[{—u + mf)/mf] and N2 = N[{—t + mf)/mf] for incoming partons i = 1,2, 
Ci = Cp = {N^ — l)/{2Nc) for quarks, with N^. the number of colors, and Cj = Ca = N^. for 
gluons. 

The second exponent resums collinear and soft gluon emission from the outgoing massless 
partons, if any (none in associated production), in the hard scattering and is given in the 
MS scheme by 
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N{s/m'f), and Bf'^ equals 3Cp/A for quarks and (3o/A for gluons, with (3q 



Here A^^ 

(llCyi—2nj)/3 the lowest-order term in the expansion of the /^-function, P{as) = fidlng/dfi = 
~/5otts/(47r) + (9(a^), where is the number of hght quark flavors. 

In the third exponent, the parton anomalous dimensions are given by 7g/g(A^i, as) = 
(a,/7r)[3C^/4 - CplnNi] + 0{al) for quarks and -fg/g{Ni,as) = (a,/7r)[/5o/4 - C^lnA/'i] + 
0{al) for gluons, and /xj? is the factorization scale. 

In the fourth exponent the (3 function is as described above, the constant da^ = or 
1 if the leading order cross section is of order q;° {t and s channels) or a] (associated tW 
production) respectively, and fiR is the renormalization scale. 

fffifi aj-g f^i^Q hard-scattering functions for the scattering of partons and fj, while S-^^'^^ 
are the soft functions describing noncollinear soft gluon emission We use the expansions 
H = + (a^/7r)if(i) + 0{al) and S = 5^ + (a,/7r)5(i) + 0{a^^). Also N = Ne^^, with 
je the Euler constant. At lowest order, the product of if-^'-^j and S^^^^ reproduces the Born 
cross section for each partonic process, = a'^'^'^ H^^^ S^^\ 

The evolution of the soft function follows from its renormalization group properties and 
is given in terms of the soft anomalous dimension ^ [211 1221 [31]. We expand Ts = 



as/'n')T^s^ + C(a^). The one-loop term T^g' is determined through the explicit calculation. 
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in the eikonal approximation, of the one-loop diagrams involving eikonal vertex corrections as 
well as top quark self-energies as shown in Figs. 4-7 (for the requisite integrals in dimensional 
regularization see plj). Here the partons in the scattering are represented by eikonal lines 
that connect in a color vertex. Explicit expressions for the soft anomalous dimensions in the 
various channels are given in the next section. 
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Figure 4: One-loop eikonal vertex corrections to the soft function for the t-channel diagram 
in single top quark production. 
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Figure 5: One-loop eikonal vertex corrections to the soft function for the s-channel diagram 
in single top quark production. 

By expanding the resummed cross section, Eq. fl2.1l) . in powers of we derive fixed- 
order corrections in the perturbative series that do not involve the prescription ambiguity 
that a fully resummed cross section entails [23]. In the n-th order corrections, the leading 
logarithms (LL) are those with / = 2n — 1 while the next-to-leading logarithms (NLL) are 
those with I = 2n — 2. In this paper we calculate next-to-leading order (NLO), next-to-next- 
to- leading order (NNLO), and NNNLO soft-gluon threshold corrections at NLL accuracy, i.e. 
at each order including both leading and next-to-leading logarithms, and also consistently 
including terms that involve the factorization and renormalization scales and ( constants, 
which arise from the inversion from moment space where the resummation is performed back 
to momentum space. We denote these corrections as NLO-NLL, NNLO-NLL, and NNNLO- 
NLL, respectively. Full details are given in the next section (see also Refs. [50]). 
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Figure 6: One- loop eikonal vertex corrections to the soft function for the associated tW 
production diagram in single top quark production. 
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(a) (b) (c) 

Figure 7: Top-quark eikonal self-energy one-loop corrections in single top quark production: 
(a) t channel; (b) s channel; (c) associated tW production. 

3 NNNLO soft-gluon corrections 



We now proceed with the calculation of the NLO, NNLO, and NNNLO soft-gluon corrections 
at NLL accuracy by expanding the resummed cross section, Eq. (12. ip . In our derivation 
of these corrections we follow the general techniques and master formulas presented in Ref. 

The NLO soft-gluon corrections at NLL accuracy for the processes are 



dt du 



TT 



C3 



ln(s4/mj) 



S4 



+ C2 
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S4 



(3.1) 



where Fb is the Born term defined in the Appendix, Eq. ( lA.Op . 

For the t and s channels the LL coefficient is c^'^ = 3Cf while for associated tW produc- 
tion it is cf^ = 2{Cf + Ca). 

The NLL coefficient, C2, can be written as C2 = T2 + C2 , where T2 represents the scale- 
independent part of C2 and C2 has all the scale dependence. For the t and s channels, 
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while for associated tW production, 
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with mw the ly-boson mass, and 



utw 

^2 



-(C^ + Ca) In 



(3.5) 



The term Rer'^^"* denotes the real gauge-independent part of the one-loop soft anomalous 
dimension (gauge-dependent terms cancel out in the cross section). A one-loop calculation 
gives for the t channel (Figs. 4, 7a) 



R^V'f = Cf 



In — + In 
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(3.6) 



for the s channel (Figs. 5, 7b) 
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m — 
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(3.7) 



and for associated tW production (Figs. 6, 7c) 



ReVf *^ = Cf In 



'mf — t\ Ca, fm'f — u' 
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2 
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2 



The coefficient Ci in Eq. (13. ip represents the scale-dependent part of the ^(54) corrections. 
For the t and s channels 



C^ln y 

and for associated tW production 
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(3.10) 



Note that we do not calculate the virtual corrections here. Our calculation of the NLO 
soft-gluon corrections includes the full leading and next-to-leading logarithms and is thus a 
NLO-NLL calculation. 
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We next calculate the NNLO soft-gluon corrections. In the t and s channels the correc- 
tions take the form 



t,s 
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(3.11) 



where C2 = vr76 and Cs = 1.2020569- 
values given previously for each channel. 

For associated tW production the corrections are 



and where for C3, C2, T2, C2, and we use the 
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(3.12) 



where for C3, C2, T2, C2, and we use the values given previously for the tW channel. We 
see that this is very similar to the expression for the t and s channels, Eq. f l3.1ip . differing 
only by some /?o terms in the [In2(s4/m^)/s4]+ and [1/54]+ terms. 

We note that only the leading and next-to-leading logarithms are complete, i.e. are 
fully known. Hence this is a NNLO-NLL calculation. Consistent with a NLL calculation 
[30] we have also kept all logarithms of the factorization and renormalization scales in 
the [In(s4/mf)/s4]+ term, and squares of logarithms of the scales in the [1/54]+ term. We 
have also kept C2 and Cs terms that arise in the calculation of the soft corrections when 
inverting from moments back to momentum space [231 ED] • This includes all C terms in the 
[In(s4/m^)/s4]+ and [1/54]+ terms, and C terms multiplying logarithms of the scales in the 
5(54) term. 

At NNNLO for all channels the corrections take the form 
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where, for the t and s channels, X3 = PQC3/I2 — CfPo/24:, X2 = (/3o/8)c3 ln(/i^/mj), 
Xf = -(/?o/4)c^ln(/iym2)-C^(/5o/4)ln2(/iym?), = (/3o/4) C2 C3 -C,.(/3o/8)C2, Xi^'^ = 
— C2C3C2, and Xq'^ = (/3o/8)C2C3 ln(/iym^). The terminology is taken from Ref. [30] • For 
associated tW production the expressions for X^, X^^, and X^'' are the same, but for 
the other variables we have X|^ = /3oC3/12, xf'^ = -{(3o/4:)c^Hl^l/m'l) - {Cp + 
Ca){Po/^) ); and Xf * = (/?o/4) C2 C3. We note that only the leading and next-to- 

leading logarithms are complete. Hence this is a NNNLO-NLL calculation. Consistent with 
a NLL calculation pH [30] we have also kept all logarithms of the factorization and renor- 
malization scales in the [In^ (54/771^ )/s4]+ term, squares and cubes of logarithms of the scales 
in the [In^(s4/mj)/s4]+ term, and cubes of logarithms of the scales in the [In(s4/m^)/s4]+ 
term. We have also kept (2 and ^3 terms that arise in the calculation of the soft corrections 
when inverting from moments back to momentum space [2ll |30]. This includes all ( terms 
in the [In^(s4/mj)/s4]+ and [In^(s4/m^)/s4]4. terms, ( terms multiplying logarithms of the 
scales in the [In(s4/m^)/s4]+ term, ( terms multiplying squared and cubed logarithms of the 
scales in the [1/54]+ term, and ( terms multiplying cubed logarithms of the scales in the 
5(54) term. 



Single top quark production via the t channel at the 
Tevatron 



We now convolute the partonic cross sections in Section 3 with parton distribution functions 
(pdf ) to obtain the hadronic cross section in pp collisions at the Tevatron. Details on the 
kinematics are given in the Appendix. We use the MRST2004 NNLO pdf [32] throughout. 
We also use standard values for the various electroweak parameters in the calculations (see 
the Appendix). 

We begin with the t-channel. The dominant processes (with percentage contribution to 
the cross section) are ub — > dt (65.7%) and db ut (21.4%). Additional processes involving 
only quarks are cb st (2.7%) and the Cabibbo-suppressed ub st (3.6%), cb dt (0.15%) 
and us dt (0.4%); the contributions from even more suppressed processes {ub bt, 
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cb — > bt, ud — > dt, etc.) are negligible. Additional processes involving antiquarks and quarks 
are sb — *• ct (4.4%) and the Cabibbo-suppressed db — >• ct (1.2%), sb — > ut (0.2%) and ds — > ut 
(0.14%); the contributions from even more suppressed processes {ss — > ct, dd ut, sd — > ct, 
etc.) are negligible. 



t channel 


LO 


NLO approx 


NNLO approx 


NNNLO approx 


mt = 170 


1.131 


1.150 


1.177 


1.193 


mt = 172 


1.091 


1.113 


1.139 


1.155 


= 175 


1.035 


1.060 


1.085 


1.100 



Table 1: The leading-order and approximate higher-order cross sections for top quark pro- 
duction in the t channel in pb for pp collisions with =1.96 TeV and = 170, 172, and 
175 GeV. We use the MRST2004 NNLO pdf and we set /i^ = /i^ = mt. 



In Table 1 we give results for the leading-order (LO) cross section and for the approximate 
NLO, NNLO, and NNNLO cross sections that include the threshold corrections at NLL 
accuracy at each order (i.e. NLO-NLL, NNLO-NLL, NNNLO-NLL). We set jjp = = mt 
and show results for three different top-quark mass values. We use the same pdf set for all 
results because we are interested in the size of the terms at each order in the perturbative 
calculation with all other things held constant. Note that the soft-gluon corrections are 
relatively small for this channel. This is also true for the exact NLO corrections ^15j. In 
fact, the approximate NLO cross section is less than 2% larger than the exact NLO result. 

The most recent value for the top quark mass from the Tevatron is mt = 171.4 ±2.1 
GeV |33], so it is interesting to calculate the cross section for this specific mass value. The 
NNNLO approximate cross section at 171.4 GeV is 1.17 pb. If we match this to the exact 
NLO cross section, then the matched cross section (i.e. exact NLO plus NNLO and NNNLO 
threshold corrections) is 

^i-channei^^^ = 171.4 ± 2.1 GeV) = 1.15i°:°? ± 0-04 ± 0.06 pb , (4.1) 

where the first uncertainty is due to the scale dependence, the second is due to the mass, 
and the third is due to the pdf. Adding these uncertainties in quadrature we find that 
^t-channei^^^ _ ]^7]^_4 -J- 2.1 GeV) = 1.15 ± 0.07 pb. A few more remarks are in order 
regarding the uncertainties. The scale uncertainty results by varying the scale between 
mt/2 and 2mt. Although this is a standard procedure, theoretically it is not unambiguous. 
The mass uncertainty in the cross section is found by using the experimentally determined 
uncertainty in the mass of the top quark (±2.1 GeV). Regarding the pdf uncertainty, since 
the MRST2004 densities do not come with errors, we use instead the pdf uncertainty from 
the MRST2001E NLO pdf [34] (the two sets give very similar values for the cross section 
and the latter set also provides pdf uncertainties). 

It is also important to provide the matched cross section for mt = 175 GeV, a top quark 
mass value that has been used widely in cross section calculations. As shown in Table 1, 
the NNNLO approximate cross section at 175 GeV is 1.100 pb. After matching to the 
exact NLO cross section we find ^^-^^^'^'^^^(mt = 175 GeV) = 1.08lo;oi ± 0.06 pb, where the 
first uncertainty is due to the scale dependence and the second is due to the pdf, and no 
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Single top at Tevatron t-channel S =1.96 TeV |-l=in^ 




uncertainty is considered for the mass. Adding these uncertainties in quadrature we find 
^t-channei(^^^ = 175 GeV) = 1.08 ± 0.06 pb. 

In Fig. [8] we plot the cross section for single top quark production at the Tevatron with 
VS = 1.96 TeV in the t channel using the MRST2004 NNLO parton densities and setting 
both the factorization and renormalization scales to a common scale /i = mj. We plot the 
LO cross section and the approximate NLO, NNLO, and NNNLO cross sections at NLL 
accuracy. 

The relative contribution from higher orders is shown in Fig. [9] where the K factors, 
defined as ratios of the higher-order cross sections to LO, are shown. The K factors in the 
t channel are small, showing that the corrections do not provide a big enhancement to the 
cross section. 

In Fig. [To] we plot the scale dependence of the cross section with rrit = 175 GeV. We set 
the factorization scale equal to the renormalization scale and vary this common scale fi over 
two orders of magnitude. We show results for the LO and higher-order cross sections. 

In general the factorization and renormalization scales are independent and it is interest- 
ing to see the dependence of the cross section separately on each scale. In Fig. [11] we plot 
the factorization scale, ^p, dependence of the cross section, while setting the renormalization 
scale fiR = rrit. We vary hf over two orders of magnitude. At LO the result is exactly the 
same as in Fig. [TD] because at LO there is no /ir dependence. However at higher orders the 
difference is noticable. 

In Fig. [12] we plot the renormalization scale, /zk, dependence of the cross section, while 
setting the factorization scale = ^t- We vary fiji over two orders of magnitude. At LO 
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Figure 9: The K factors for single top quark production at the Tevatron in the t channel. 
Here fi = fip = fiR = mt- 
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Figure 10: Scale dependence of the single top quark cross section at the Tevatron in the t 
channel. Here fi = fip = ^r. 
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Figure 11: Factorization scale dependence of the single top quark cross section at the Teva- 
tron in the t channel. Here the renormalization scale is held fixed, = rrit = 175 GeV. 



the result is exactly fiat because, as noted above, at LO there is no fiR dependence. At higher 
orders there is a fiR dependence, though it is fairly fiat compared to the fip dependence in 

Fig. m 



5 Single top quark production via the s channel at the 
Tevatron 

We continue with the s channel. The dominant process (with percentage contribution to the 
cross section) is ud —>■ ht (97.4%). Additional processes are cs — > ht (1.1%) and the Cabibbo- 
suppressed us — > ht (1.1%), cd ht (0.26%), and ud — > st (0.17%); the contributions from 
even more suppressed processes {ud — >• dt, cs — > dt, etc.) are negligible. 



s channel 


LO 


NLO approx 


NNLO approx 


NNNLO approx 


mt = 170 


0.353 


0.510 


0.555 


0.576 


mt = 172 


0.335 


0.484 


0.528 


0.547 


mt = 175 


0.310 


0.448 


0.488 


0.506 



Table 2: The leading-order and approximate higher-order cross sections for top quark pro- 
duction in the s channel in pb for pp collisions with = 1.96 TeV and mt = 170, 172, and 
175 GeV. We use the MRST2004 NNLO pdf and we set ^f = Ij^r = mt. 



In Table 2 we give results for the LO cross section and for the approximate NLO, NNLO, 
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Figure 12: Renormalization scale dependence of the single top quark cross section at the 
Tevatron in the t channel. Here the factorization scale is held fixed, hf = rrit = 175 GeV. 



and NNNLO cross sections with the threshold corrections at NLL accuracy. The soft-gluon 
corrections are relatively large for this channel, in stark contrast with the results we found in 
the t channel. This is also true for the exact NLO cross section [15]. That the behavior of the 
two channels is quite different should not be surprising: the kinematics and the color flows 
are quite different. The t channel resembles deep inelastic scattering while the s channel 
resembles the Drell-Yan process. As we will see, there are differences between the channels 
not only in the K factors but also in the scale dependence. We also note that the approximate 
NLO cross section in the s channel is only 3% larger than the exact NLO result, showing 
that the threshold corrections are dominant and provide the bulk of the QCD corrections 
and, thus, that the threshold approximation works very well. 

Again, it is interesting to provide the cross section for the specific value of the new mass 
from the Tevatron, rrit = 171.4 ±2.1 GeV. The NNNLO approximate cross section at 171.4 
GeV is 0.555 pb. If we match this to the exact NLO cross section, then the matched cross 
section (i.e. exact NLO plus NNLO and NNNLO threshold corrections) is 



channel / 



[mt 



171.4 ± 2.1 GeV) = 0.54 ± 0.02 ± 0.03 ± 0.01 pb . 



(5.1) 



where the first uncertainty is due to scale variation between mj/2 and 2mj, the second 
is due to the mass (±2.1 GeV), and the third is the pdf uncertainty, as discussed in the 
previous section. Adding these uncertainties in quadrature we find that cr^~^^^^'^'^^{mt = 
171.4 ± 2.1 GeV) = 0.54 ± 0.04 pb. 

As before, it is also important to provide the matched cross section for rrit = 175 GeV. 
As shown in Table 2, the NNNLO approximate cross section at 175 GeV is 0.506 pb. After 
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Figure 13: The cross section for single top quark production at the Tevatron in the s channel. 
Here fi = fip = fiji = rrit. 



matching to the exact NLO cross section we find o-^-^i^^'^°<=i(mt = 175 GeV) = 0.49±0.02±0.01 
pb, where the first uncertainty is due to the scale dependence and the second is due to the 
pdf, and no uncertainty is considered for the mass. Adding these uncertainties in quadrature 
we find a'-''^^'"^''\mt = 175 GeV) = 0.49 ± 0.02 pb. 

In Fig. [13] we plot the cross section for single top quark production at the Tevatron with 
=1.96 TeV in the s channel using the MRST2004 parton densities and setting both the 
factorization and renormalization scales to /i = m^. We plot the LO cross section and the 
approximate NLO, NNLO, and NNNLO cross sections at NLL accuracy. 

The K factors are shown in Fig. O They are quite large, thus showing that the 
corrections provide a big enhancement to the cross section. This behavior is to be contrasted 
with the t channel where the corrections are small. We also note that the i^-factors in the 
s channel are fairly constant over the top-quark mass range shown. As seen in the plot, 
the NLO corrections provide a 45% increase of the LO cross section, the NNLO corrections 
provide an additional 12%, and the NNNLO corrections a further 6%. 

In Fig. [15] we plot the scale dependence of the cross section with rrit = 175 GeV. We set 
the factorization scale equal to the renormalization scale and vary this common scale /i over 
two orders of magnitude. We show results for the LO and higher-order cross sections. 

Again, we note that the factorization and renormalization scales are theoretically inde- 
pendent, and it is interesting to see the dependence of the cross section separately on each 
scale. In Fig. [T^we plot the factorization scale, fip, dependence of the cross section, while 
setting the renormalization scale fiR = rrit. We vary fip over two orders of magnitude. At 
LO the result is exactly the same as in Fig. [15] because of the lack of /i/j dependence. 
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Figure 14: The K factors for single top quark production at the Tevatron in the s channel. 
Here ^ = ftp = fiR = rrit. 
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Figure 15: Scale dependence of the single top quark cross section at the Tevatron in the s 
channel. Here ^ = iip = I^r- 
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Figure 16: Factorization scale dependence of the single top quark cross section at the Teva- 
tron in the s channel. Here the renormalization scale is held fixed, /ir = rrit = 175 GeV. 

In Fig. [T7| we plot the renormalization scale, fin, dependence of the cross section, while 
setting the factorization scale fip = nit. We vary fiji over two orders of magnitude. At LO 
the result is exactly flat, as noted before. At higher orders there is a fin dependence and it 
is less flat than for the t channel. 

6 Associated tW production at the Tevatron 

Associated tW production proceeds via bg tW~ (Cabibbo-suppressed contributions from 
sg — >■ tW~ and dg —>■ tW~ are negligible). The kinematics and color flow for this process 
are rather different from both the t and s channels since we have a real W boson in the flnal 
state and a quark-gluon vertex at lowest order. 

The LO cross section is 0.070 pb for = 175 GeV, which is a rather small cross section 
and thus makes this channel relatively unimportant at the Tevatron even after we have added 
the threshold corrections. 

In Table 3 we give results for the LO and the approximate NLO, NNLO, and NNNLO 
cross sections that include the threshold corrections at NLL accuracy at each order. The 
soft-gluon corrections are relatively large for this channel, even more than in the s channel. 

Again, it is interesting to give a result for the speciflc value of the new mass from the 
Tevatron, = 171.4 ±2.1 GeV. The NNNLO approximate cross section at 171.4 GeV is 



where the flrst uncertainty is due to scale variation between mt/2 and 2mt, the second is due 



a^^imt = 171.4 ± 2.1 GeV) = 0.14 ± 0.02; 



,+0.01 
-0.02 



±0.02 pb. 



(6.1) 



17 
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Figure 17: Renormalization scale dependence of the single top quark cross section at the 
Tevatron in the s channel. Here the factorization scale is held fixed, hf = rrit = 175 GeV. 



tW production 


LO 


NLO approx 


NNLO approx 


NNNLO approx 


mt = 170 


0.080 


0.122 


0.140 


0.146 


mt = 172 


0.076 


0.116 


0.133 


0.139 


= 175 


0.070 


0.107 


0.122 


0.127 



Table 3: The leading-order and approximate higher-order cross sections for associated tW 
production in pb for pp collisions with = 1.96 TeV and mt = 170, 172, and 175 GeV. 
We use the MRST2004 NNLO pdf and we set = = mt. 



to the mass (±2.1 GeV), and the third is the pdf uncertainty, as discussed previously. Adding 
these uncertainties in quadrature we find that a'^{mt = 171.4 ± 2.1 GeV) = 0.14 ± 0.03 pb. 
We note that since an exact NLO result is not available we do not provide a matched cross 
section here. 

As before, it is also important to provide the cross section for mi = 175 GeV. As shown in 
Table 3, the NNNLO approximate cross section at 175 GeV is 0.127 pb. Including uncertain- 
ties we have (T*^(mj = 175 GeV) = 0.13 ±0.02 ±0.02 pb, where the first uncertainty is due to 
the scale dependence and the second is due to the pdf, and no uncertainty is considered for 
the mass. Adding these uncertainties in quadrature we find a'^imt = 175 GeV) = 0.13±0.03 
pb. 

In Fig. [18] we plot the cross section for single top production at the Tevatron with 
= 1.96 TeV via associated tW production using the MRST2004 parton densities and 
setting both the factorization and renormalization scales to fi = mt. We plot the LO cross 
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Figure 18: The cross section for associated tW production at the Tevatron. Here jj, = jip = 
fJ'R = "bi- 



section and the approximate NLO, NNLO, and NNNLO cross sections at NLL accuracy. 

The relative contribution from higher orders is shown in Fig. [191 The K factors are 
quite large, showing that the soft-gluon corrections provide a big enhancement to the cross 
section, and fairly constant over the top-quark mass range shown. As seen from the plot, 
the NLO corrections provide a 53% increase of the LO cross section, the NNLO corrections 
provide an additional 22%, and the NNNLO corrections a further 8%. However, the overall 
cross section still remains rather small and we will not study it further here. 

7 Conclusion 

We have studied single top quark production at the Tevatron via all Standard Model partonic 
processes. We have presented the resummation of the threshold soft-gluon corrections and 
their expansion to provide higher-order approximate cross sections through NNNLO. 

The results differ a lot among channels. In the t channel the soft-gluon corrections 
are small and hence do not greatly change the LO result. Including soft gluon corrections 
through NNNLO at NLL accuracy and matching with the exact NLO cross section, our best 
estimate of the cross section is o-*-channeij^^^ _ iji ji^ j- 2.1 GeV) = 1.15 ± 0.07 pb, where the 
uncertainty indicated includes the scale dependence, the uncertainty in the top quark mass, 
and the pdf uncertainty. 

In the s channel, however, the corrections are significant. The NLO approximate cross 
section is an excellent approximation to the full NLO result, thus showing that the thresh- 
old approximation works very well, while the higher-order contributions provide further 
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Figure 19: The K factors for associated tW production at the Tevatron. Here jj, = (ip = 
fiR = rrit. 



171.4 ± 2.1 GeV) = 0.54 ± 0.04 pb. 

The renormahzation and factorization scale dependence of the cross section for both t 
and s channels were displayed both separately and also setting the scales equal to each other. 

The threshold corrections for associated tW production are large as a percentage of the 
LO cross section (large K factors), however even with those the total cross section is still 
rather small for this channel. We find a*^(mt = 171.4 ± 2.1 GeV) = 0.14 ± 0.03 pb. 

Finally, we note that the corresponding cross sections for single antitop quark production 
at the Tevatron are identical to those for single top quark production as presented in this 
paper. 
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enhancements. Our best estimate for the cross section in this channel is o" 



s— channel 
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Appendix 



Here we provide some details on the kinematics of the calculation. We consider the process 
Pi + P2 ^ Ps + PA:, where the particles have been represented by their momenta pi. If pi 
and p2 refer to incoming partons (quarks or gluons) then rrii = and 1712 = 0. The partonic 
kinematical invariants are s = (pi+P2)^, t = (pi— Ps)^, u = {P2—P3Y1 = s+t+u—m'^—ml. 
The hadronic kinematical invariants are S = {phi +Ph2Y, T = {phi —PzY^ U = {ph2 —PsY, 
5*4 = S + T + U — m\ — m\, where hi and /12 are the hadrons (protons and antiprotons at the 
Tevatron) corresponding to partons pi and p2] thus, pi = XiPhi and p2 = X2Ph2- The relations 
between partonic and hadronic quantities are as follows: s = X1X2S, t — m\ = Xi{T — mg), 
u — m| = X2{U — rri^). Then 

-4 = --{l-xiy ^-{1-X2)^ ^ + 0((1-Xl)(l-X2)). (A.l) 

b s s s 

If we denote by a the hadron-level cross section and by a the parton-level cross section, 
then we have the relations 

a = J dxidx2(f){xi) (f){x2) & 
da r , dt du da 



dTdU 

where (/>(a;) are the parton distributions. Now dt = xidT and du = X2dU, so 

'dXidX2XiX2<^iXi)<^ix2) 



dT dU J dt du 

dTdUdx2ds4 — ^—7^^ 2 <^(^2) ^r-^- , (A. 3) 



where we used the relation 



_ SA-ml + ml-X2{U -ml) 
X2S + T — ml 

Writing explicitly the kinematical limits on the integrals, we get the precise expression 
for the hadronic cross section 

'^VhiPh2^PsvAS) = dT dU dx2 ds^ 



X ^ l'^' 2 <Pi^^) <Pi^2) (A.5) 



where Xi is given in Eq. (lA.4p . and the hmits on the integrals are 



T^: = -l{S-ml-ml)± y{S-ml-miy-4mlml, (A.6) 



mis 

h 
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Umax = ml+ ^'^ ^ , Umin = -S -T + ml + mi, (A. 7) 

J TTlo 



and s^rnax = X2{S + U - m^) + T - ml. 

The Born-level differential partonic cross section is written as 

dtdu - ^v.v.^v.vA^^) - i6vrs2 ^^'^^ ' ^^'^^ 

where the amplitude squared has to be specified for the partonic process studied. 
For the t-channel process uh —>■ dt we have 



Ati'^o^VHV^^^ s(s -ni- 



2^ 



\Mu,^,,\' = ^-T47^"%r ' % . (A.10) 
sm tlw [t — m^rY 



For the t-channel process dh ut we have 



sin^ 6w{t — ml^Y 
For the s-channel process ud bt we have 



in^a^lV^, t{t - m 



For the associated production process bg tW~ we have 

\M, ? = f ^ _ + , (A.13) 

bg^tw Sm^,/ sin^ 6'vK \Mi wis / ' 

where A 1 = -^2(5 -m^ - m^)(m^ + m^/2) - (ti/2)(-2m^ + m^mj +m^) -2u2m^(2m^ + 
m^), = -ti(-m^ + m^)ml, - {u2/2)tim^ - Ui{u2/2)m^ - sm^^ml - (s/2)m^, A3 = 
— Mi(s/4)(2m^ + mj) where a = e^/(47r) and is the strong coupling. 

The expressions for the corresponding processes involving charm quarks and Cabibbo- 
suppressed contributions are similar. 

The electroweak parameters used in the above expressions and throughout the paper are 
the following [33]: 

mw = 80.403 GeV, 

Gp = = 1.16637 X 10-5 GeV-2, 

sm^Owlz = 0.23122, 
e = gwsmOw = 0.31404, 



a = eV(47r) = 7.848 x 10 



-3 



For the CKM matrix elements we have used [35] 
Vud = 0.97383, Vus = 0.2272, Kb = 0.00396, 
= 0.2271, Ks = 0.97296, Kb = 0.04221, 
Vtd = 0.00814, Vts = 0.04161, K = 0.9991. 
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